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(N ' Abstract 

Recently Tang and Boozer [Phys. Plasmas (2000)], have investigated the anisotropics in 



43 

c3 



magnetic field dynamo evolution, from local Lyapunov exponents, giving rise to a metric tensor, 
in the Alfven twist in magnetic flux tubes (MFTs). Thiffeault and Boozer [Chaos(2001)] 
have investigated the how the vanishing of Riemann curvature constrained the Lyapunov 
exponential stretching of chaotic flows. In this paper, Tang-Boozer-Thiffeault differential 

m ! 

geometric framework is used to investigate effects of twisted magnetic flux tube filled with 

cn ■ 

helical chaotic flows on the Riemann curvature tensor. When Frenet torsion is positive, the 
Riemann curvature is unstable, while the negative torsion induces an stability when time 

O ' 

t — > oo. This enhances the dynamo action inside the MFTs. The Riemann metric, depends 
on the radial random flows along the poloidal and toroidal directions. The Anosov flows 
has been applied by Arnold, Zeldovich, Ruzmaikin and Sokoloff [JETP (1982)] to build a 
uniformly stretched dynamo flow solution, based on Arnold's Cat Map. It is easy to show 
that when the random radial flow vanishes, the magnetic field vanishes, since the exponential 
Lyapunov stretches vanishes. This is an example of the application of the Vishik's anti-fast 
dynamo theorem in the magnetic flux tubes. Geodesic flows of both Arnold and twisted MFT 
dynamos are investigated. It is shown that a constant random radial flow can be obtained 
from the geodesic equation. Throughout the paper one assumes, the reasonable plasma 
astrophysical hypothesis of the weak torsion. Pseudo-Anosov dynamo flows and maps have 
also been addressed by Gilbert [Proc Roy Soc A London (1993)]. 

PACS numbers:02.40.Hw:differential geometries. 91. 25. Cw-dynamo theories. 
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I Introduction 



To talk about the importance of the developement of the concepts of Riemannian geom- 
etry [1] to physics, would be enough to mention the advances in Einstein's gravitational 
theory, more well known as the general theory of relativity [2]. More recently, uses in 
other parts of physics such as string theory of defects and solid state theory [3] have also 
experienced the same success. Yet more recently, Thiffeault , Tang and Boozer, on a series 
of papers [4, 5, 6, 7] have investigate the constraints experienced by chaotic flows, when 
the Riemann curvature tensor vanishes. In this paper, a slightly more general application 
of these ideas are applied to helical chaotic flows inside the twisted magnetic flux tubes 
[8, 9] is given. The role of random radial flows in Lyapunov exponential stretching, so 
fundamental to the kinematic dynamo problem, is investigated throughout the paper. 
The investigation of Anosov diffeomorpisms [9] , an important mathematical tool from the 
theory of dynamical systems, has been often used, in connection with the investigation 
of dynamo flows and maps [9] such as the Arnold's Cat Map [9] on the torus, useful in 
mixing [10] problems in the physics of fluids [4]. One of the main properties of the Anosov 
maps is that they yield Lyapunov exponential of the chaotic exponential stretching, which 
are constant everywhere [5]. In their paper, Thiffeault and Boozer [5], obtained geomet- 
rical constraints in two and three dimensions, with coordinates x e R 2 and x' G R 3 . 
Consequences of the vanishing of the Riemann curvature tensors, or Riemann-flat condi- 
tions in the mathematicians jargon, includes constraining the modification of the growth 
rate of magnetic fields in the kinematic dynamo problem or the mixing properties of 
chaotic advection-diffusive equation. They conclude that, the Lyapunov exponents are 
locally small in regions where, Riemann curvature of the stable manifold, is large. Though 
Euclidean metrics were used in their work to allow for the vanishing of the Riemann cur- 
vature, here this hypothesis is relaxed, which still allows one to obtain similar results for 
the flux tube endowed with chaotic helical flows. The Riemann and Ricci curvatures are 
computed in terms of the Lyapunov exponential stretching and establishing the conditions 
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on the Frenet torsion of the tube axis to obtain positive Lyapunov exponents which leads 
to the fast dynamo action for chaotic flows permeatting the tube. 

This last condition is fundamental for dynamo action [8]. Random radial flows [10, 11] are 
shown to be fundamental, for the chaotic flows to experience a non- vanishing exponential 
stretching. Since these MFTs are fundamental in astrophysical plasmas and solar chaotic 
plasmas in tokamaks and stellarators, as well in the Perm torus of liquid sodium dynamo 
[12], the better understanding of chaotic flow, seems to be useful to future developement 
of plasma physics and dynamo theory. Therefore a good understanding of the behaviour 
of these rotating torus flows is of utmost importance, in future dynamo experiments. 
This paper is organised as follows. In section 2, special types of chaotic flows in twisted 
MFTs are presented such as the torus automorphisms, called Arnold's Cat Map, which 
led Arnold [9] to build a uniform stretching dynamo solution of self-induction equation. 
In this case though the Jacobian differs from the Arnold, Anosov map one, the same 
eigenvalues are obtained, and therefore the dynamo chaotic flow is stretched along one 
direction and compressed along the other. Section 3 generalizes the Anosov flows to the 
case where the stretch along toroidal direction is not constant and therefore the Lyapunov 
metric of the tube leads to the computation of Riemann and Ricci curvature in the twisted 
MFTs metric. This kind of twisted geometry is then used to established constraints along 
the section to established the stability relation to the torsion of the tube, which in turn 
is proportional to the twist. Section 4 presents the discussion and conclusions. 
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II Hyperbolic dynamo maps and flows in MFTs 



This section contains a brief review and some new material on hyperbolic flows and 
dynamo maps. While is not the main section of the paper, it is supposed to pave the 
way to the next section where the main results are given. Let us start by defining what 
exactly is an hyperbolic map. Anosov maps are actually an example of hyperbolic map 
with extra constraints which one shall not bother here. An hyperbolic map on a torus 
may be defined as an automorphism where the determinant of its Jacobian is equal to one. 
This constraint shall be used in this section, in analogy with Thiffeault and Boozer [4] 
constraints to place constraints on the dynamos and Riemannian geometry of the tubes. 
Let us consider the following two-dimensional twist map [13] Jacobian matrix 

j t wist-(; 2) (iu) 

Here K is the toroidal stretching along the tube metric [14] given by 

dl 2 = dr 2 + r 2 d9 R 2 + K 2 (r,s)ds 2 (11.2) 

where of course here, the internal cross- section radius r is constant for the two dimensional 
Jacobian and the twist transformation is 

6(s) :=e R - J r(s)ds (11.3) 

where the integral in this expression is the total torsion. Here, since one is considering 
helical chaotic flows, the Frenet torsion and curvature of the tube axis are constants and 
fixed to be equal. As one shall note in the next section the introduction of coordinate r 
into the stretching K(r, s) = l — rK(s)cos6 allows one to has Riemann and Ricci curvature 
that do not vanish identically. Of course here, one have considered that K is a constant 
or uniform stretching so important in chaos [14]. A small lemma can be proved here by 
constrained the Jacobian transformation to be hyperbolic. Thus 

lemma: Hyperbolic maps with constant stretch in the twist MFTs, leads necessarilly to a 
thin tube. Proof: This comes immeadiatly from the determinant of matrix (II. 1) equals 
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to one 

^[ J twist] = K = 1 (II.4) 

which by examining the expression for K(r, s) shows that K(0, s) = 1 is the thin tube 
approximation, as one wishes to prove. Let us now compute the eigenvalue problem. Just 
to show that the twist tube Jacobian above is a general case of the torus automorphism by 
Arnold, let us consider a special transformation where the twist is fixed while the radius 
is contracted to half his size and the stretch is double, similar to what happens in the 
stretch-twist-fold (STF) dynamo mechanism proposed by Vainshtein and Zeldovich [15]. 
Mathematically this transformation is 

, 1 



2 

s' = 2s 



(II.5) 
(II.6) 



while & = Or. Besides the torsion To is normalized to one, since one does not want to 
deal with untwisted tubes. Therefore the three-dimensional Jacobian 

(K x \ 



'3D = 1 

V K J 

where K\ is the stretching (contraction) of the radial direction, takes the form 

(\ \ 



(II.7) 



'3D 



(II.8) 



1-1 
\0 2 / 

By taking the determinant of this Jacobian one obtains it is equal to one, or hyperbolic, 
and besides the eigenvalue problem 



Det[J 3B - XI] = 



(II.9) 



where 1 is the unity matrix 





n 





o\ 


1 = 





1 
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(11.10) 
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A being the eigenvalue. The determinant leads to the following second-order algebraic 
equation 

A 2 -3A + 1 = (11.11) 
Note that this equation leads to the Arnold eigenvalues 

a ± = ^ cn.12) 

where A+ > represents the stretching, while A_ < is the contraction of the flow. In the 
Arnold's case [13] the stretching is uniform. In the next section, one shall be concerned 
with the non-uniform stretching where the stretching Riemann metric factors, do depend 
upon the coordinates. Another interesting lemma can be proved for the three-dimensional 
matrix, as 
lemma: 

The hyperbolicity condition applied to twist map in R 3 , given by the above Jacobian, 
implies that the magnetic flux tubes stretching and contraction obeys the following con- 
ditions: Kq = K\~ x , which tell us that the stretching is the inverse of contraction as 
foreseen by our physical intuition. 
Proof: 

The proof is trivial, and obtained by applying the DetJ = 1 condition to the three- 
dimensional Jacobian Jqj) above. This yields that KqKi = 1 as we wish to prove. 
The domain of the flow is given by a compact Riemannian manifold represented as the 
product of a torus T 2 and the closed interval [0, 1] of < z < 1 [9]. This results in the 
Arnold Riemann metric 

ds 2 = e~ Xz dp 2 + e Xz dq 2 + dz 2 (II. 13) 

which represents the stretching and contraction in distinct Euclidean directions p and 
q , respectively induced by the eigenvalues A± = 3± /^ also corresponding to magnetic 
eigenvectors. The stationary dynamo flow considered by Arnold was given by the simple 
uniformly stretching flow v = (0, 0, v) in (p, q, z) coordinates. More recent attempts to 
build a dynamo action by making use of compact Riemannian geometry includes the case 
of the fast dynamo of Chiconne and Latushkin [16], and the conformally stretched fast 
dynamo by Garcia de Andrade [17]. Anosov map in Arnold's cat map, is given by the cat 
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dynamo Jacobian matrix [4] 

J cat= ^ )] ( IL14 ) 

Note that this Jacobian is distinct to the previous Jacobians but satisfies the hyperbolicity 
condition and leads to the same eigenvalues as before. Let us now compute the Riemann 
curvature of the Arnold metric 

#1212 = A 2 (11.15) 

i?i3i3 = -A 2 (11.16) 

#2323 = "A 2 (11.17) 

for the Riemann curvature tensor, while the Ricci tensor components are 

r 33 = r = 2A 2 (11.18) 

the determinant of the metric g — 1, implies that the flow is incompressible as considered 
previously by Arnold. Note as well, that, as usual in Anosov systems the Riemann 
curvature is constant. 

Ill Riemann curvature from Lyapunov exponents and 
dynamo action 



In this section it is shown that in regions of weak torsion the Riemann and Ricci curva- 
tures acquire simple forms, in terms of of exponential stretching and torsion. Negative 
torsion in helical chaotic flows, leads to singularities in curvatures in the t — > oo time 
limit. A simple examination of the Jacobians considered in the last section allows us to 
see that the Riemann metric of the flux tubes vanish for a constant stretch-twist and con- 
traction. However since curvature or folding is a fundamental process in the STF dynamo 
mechanism, in this section it is shown that the consideration of a more general Jacobian 
makes the manifold acquire a curved Riemannian metric distinct from the Euclidean norm 
metrics considered so far. So, for a non-uniform stretching [14] usually found in general 
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chaotic flows the Riemann and Ricci curvatures do not vanish identically. Let us start 
computing the curvatures corresponding to the general Riemann metric of the twisted 
MFT above. With the aid of an adaptation of the tensor package the computation of the 
curvatures is neither tedious nor long and results in 

1 T W(fl) 

^ 1313 " 4(1-T rcos6) (IIU9) 

R^s = -V 0W ^ (HI.20) 
3 3 8 (l-T rcos6) y ! 

^ 2323 - 4(1-T rcos8) (HL21) 
for the Riemann curvature tensor, while the Ricci tensor components are 

4(1 - 2T rcos6 + T 2 rcos6) v ; 

1 TQ 2 rsin29 
12 ~ ~8 (1 - 2r rcos6 + T 2 rcos6) ^ ' ' 

R33 = - ~ m (HI.24) 

4 (1 - r rcos9) 

R - 1 T 2 rsin26 
22 8 (1 - 2r rcos9 + T 2 rcos6) 1 ' ; 

Most of these curvatures can be easily simplified in the weak torsion case. Let us now 
give a more dynamical character to these expressions by considering the Thiffeault-Boozer 
Riemann metric relation with the Lyapunov exponents of exponential stretching by ex- 
pressing this metric in terms of the directions t, e r and along the curved tube as 

gij = A r e r e r + A e e e ee + A s tt (111.26) 

where = r, 9, s) and Aj are the Lyapunov numbers which are all positive or null. The 
Lyapunov exponents are given by 

A~ = Um^i 1 ^) (111.27) 

The infinite symbol over the Lyapunov exponent indicates that this is a true Lyapunov 
exponent obtained as the limit of the finite-time Lyapunov exponent A«. Let us now 
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compute the values of the Lyapunov exponents which are fundamental for stretching and 
dynamos, in terms of a random radial flow 

In A • 

Xi = Jim^oof-^) (IIL28) 

<r>= J <v r >dt (111.29) 

where here one shall consider that on a finite-time period, the random flow can be con- 
sidered as approximately constant which reduces the last expression to 

< r >=< v r > t (111.30) 

Let us now compute the Lyapunov exponents of the curved Riemannian flux tubes in 
terms of the random flows as 

Inl 

\ r = lim t _ 00 (—) = Q (111.31) 
X e = limt^C-) =< v r > (111.32) 

and finally 

,lnK(r,s) —TorsinO, 
\ s = hm^oo( — ) w l«m ( ^oo( — ) ~ -r < v r > sinO (111.33) 

By following the computations by Tang and Boozer [5] on the magnetic fields by the self- 
induction equation one can say that the corresponding magnetic field components maybe 
written in terms of the Lyapunov exponents as 

B e ^ e 2Xet = e <Vr>t (111.34) 

B s w e Xst = e -ro<vr>sinet (111.35) 

while the radial magnetic field does not depend on time , which allows us to say that 
this is due to the confinement of the radial flow on the tube, actually in the force-free 
case considered by Ricca [18] the flux tube radial magnetic component B r is assumed 
to vanish. From the last two expressions it is easy to observe that the magnetic field is 
stationary in the absence of random radial flows, and there is no fast dynamo action as 
well as no stretching due to the vanishing of the Lyapunov exponents. This phenomenon 
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is actually the Vishik's anti-fast dynamo theorem [19] where no fast-dynamo action can 
be obtained, in non-stretching flows. Fast dynamo is obtained when the random flow has 
a positive average velocity. Note that in this case at least the B s toroidal field grows 
in time, while the poloidal field is spatially periodic and may even grow in time in the 
case torsion is negative. To simplify matters let us now compute the Riemann and Ricci 
curvatures for the chaotic flow tube metric 

ds 2 = dr 2 + e <Vr>t d9 2 + e~ T0<Vr>coset ds 2 (111.36) 

Note that for this expression a positive torsion shows that one of the Ljapunov exponents 
in this Riemann metric is positive (stretching) while the other is negative, representing 
contraction. The only non-vanishing component for the Riemann tensor is 

R 2323 = -^ e 2<Vr>To{1 ~ cose)t sin 2 9 (111.37) 

which shows that the Riemann curvature is unstable in the infinite time limit, if the torsion 
is positive, while it is stable if the torsion is negative. The Ricci tensor components are 

R 33 = e (i-r cos9)<v r >t R22 (111.38) 
R 22 = lcos9e- {1 - T0<Vr>t) 

The Ricci tensor is 



: 22 = -cos9e-^- T0<Vr>t} (111.39) 



R = ^ e - Xe(1 - T " )t stn 2 9 (111.40) 

In all these computations one has assumed the weak torsion approximation, which is 
very reasoble for example, in astrophysical plasmas. Note from the Riemann curvature 
expression that the Riemann space is flat when the random radial flow vanishes, and when 
the torsion To > 1 the scalar curvature R is singular, or unstable, when t — > oo. However, 
since one assumes here that the weak torsion approximation, the curvature computations 
are ingeneral unstable. Since the twist of plasma MFTs is proportional to torsion of the 
MFT axis, this is a very reasonable approximation since the twist in kink solar loops 
for examples is very weak of the order of Tw xs 10 -10 cra -1 . When torsion vanishes the 
stationary Riemann curvature is periodic and reduces to 

i? 23 23 = -\sin 2 9 (111.41) 
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Thus when the radial random flow vanishes the Riemann curvature can be spatially pe- 
riodically. Let us now compute the determinant of the Riemann metric 



This might be g = 1 in the case of incompressible chaotic flow as V.v = 0, which allows 
us to say that the chaotic helical flow inside the tube has to be compressible. Tang 
and Boozer [5], have also examine the problem of compressible chaotic flows. As a final 
observation one notes that the Frenet curvature of material lines inside the tube, increases 
as the Ljapunov exponents or metric decreases. This can be seen by a simple examination 
of the metric factor g ss = (1 — K(s)rcosO) 2 since the Frenet curvature k(s) growth implies 
that this stretching factor decreases, of course along the Ljapunov exponents since they 
depend upon the Riemann metric. 

IV Geodesic flows in chaotic dynamos 

Anosov [14] demonstrated that hyperbolic systems which are geodesic flows [20], are 
necessarilly Anosov. In this section one shows that as a consequence of the geodesic flow 
condition imposed the MFT twisted flows are non-Anosov since their Lyapunov exponents 
are not all constants. In this section two given, the first is the Anosov geodesic flow of 
Arnold dynamo, while the second addresses the example of geodesic flows on the twisted 
MFTs above. In the first case, the Arnold metric of previous section, yields the Riemann- 
Christoffel symbols as 



g = Det(gij) = r 2 (l - r rcos9) 2 



(111.42) 




9 ll (9ij,k + 9ik,j - 9jk,i) 



(IV.43) 



whose components are 




-A 




(IV.44) 




-2\z 



(IV.45) 




(IV.46) 



Substitution of these symbols into the geodesic equation 



— + r jk v j v k = o 



(IV.47) 
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yields the following equations for the geodesic flow v as 

+ r^V = o (iv.48) 



dt 

which yields 



^ _ At; V = o (IV.49) 
at 

Proceeding the same way with the remaining equations yields 

^! + A V V = (IV.50) 

and 

^! _ Afe- 2 ^^ 1 ) 2 - e 2Xz (v 2 ) 2 ] = (IV.51) 
at 

By assuming that the v 3 = constant = v 3 one obtains 

v 1 = v o e ^+v 3 ot] (IV.52) 

Since v 1 — v p — ^ and v 2 = v q = ^, one is able to integrate the geodesic equations to 
obtain the flow topology as 

p 2 + q 2 = R 2 e XLt sinh\z (1V.54) 

where := Vq\ is the Lyapunov exponent for the Arnold metric, and z := v$t. This 
geometry represents a circle expanding on time and exponentially stretched by the action 
of the Lyapunov exponent. Note that the dynamo flow velocity here, is distinct from the 
Arnold's one, since the Arnold's dynamo flow does not possesses components v p and v q 
as in the example considered in this section. Let us now proceed computing the same 
geodesic flow in the case of twisted MFT metric. Their Christoffel symbols are 

„-<?V>(l-To)t 

r 2 33 = sine (IV.55) 

e -<V r >(l-T t 

r 3 23 = sinO (IV.56) 

From these expressions one obtains the following geodesic equations 

£ = ° 
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which implies v r = v 1 = v = constant. This result is indeed important since, it shows 
that actually our hypothesis that the random radial flow is constant is feasible, and can be 
derived from a geodesic flow. To simplify the remaining equations one chooses a torsion 
To — \ , which yields 

dv 2 

— - sin6(v z f = (IV.58) 

and 

^! _ e M«v>V v 3 = (IV.59) 
dt y ' 

These two last equations together yield 

dv 3 i^, ^v 2 



, , e-3<^>*- (IV.60) 
dv 1 v A 

Some algebra yields the following result 

v _± = e (i+cose)<v r >t (IV>61) 

Thus the toroidal flow is proportional to the random flow and to the Lyapunov exponential 
stretching. Note that at t — there is an equipartition between the poloidal and toroidal 
flows while as time evolves the Lyapunov exponential stretching makes the toroidal flow 
velocity to increase without bounds with respect to the poloidal flow. This is actually 
analogous to the dynamo action. Equation (IV. 57) also shows that one the Lyapunov 
exponents is constant while the other is not and therefore the flow is a non-Anosov flow. 
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V Conclusions 



In this paper the importance of investigation the twisted torus dynamo hyperbolic chaotic 
maps and flows, has been stressed. The ideas developed here followed the path of realistic 
dynamo maps obtained from the Gilbert's [21] shearing generalization of Arnold's cat map. 
Actually the Tang-Boozer's [5] idea of pilling up the Alfven twist of magnetic flux tubes 
is fully developed geometrically by applying their Lyapunov like metric and the Riemann 
geometrical Ricca's model for the MFTs. Though Ricca has investigated only force-free 
MFTs, in this paper the fast growth of magnetic fields in MFTs, induces a dynamo 
action through the Lyapunov fast dynamo exponential stretching. It is easy to show that 
from the self-induced equation that the twisted MFTs, with zero Lyapunov exponential 
stretching leads to non-fast dynamo action in the tubes, a result which is nothing but the 
statement of anti-dynamo theorem of Vishik's [22] , which has been recently discussed in 
the context of the MFTs [23]. Computation of geodesic flows in Arnold's fast kinematic 
dynamo and in twisted MFT, shows that in both case the flow depends on the exponential 
stretching given by the Lyapunov numbers. It is shown that at the beginning of time, 
there is an equipartition between chaotic toroidal and poloidal flows, while as time evolves 
an instability occurs due to the action of Lyapunov exponential stretching. This result 
was obtained from the Euler equations by Friedlander and Vishik [24]. 
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